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I. INTRODUCTION 

While the electroweak contribution to the anomalous 
magnetic moment of the muon a M = ^(g^ — 2) can be 
calculated with sub-ppm precision, the hadronic contri- 
bution a^ ad cannot be evaluated by means of perturbative 
techniques within Quantum Chromodynamics (QCD). 
However, it is possible to relate ajj ad via a dispersion 
relation to hadronic cross section er(e + e~ — > hadrons) 
data [1]. Thus, hadronic cross section measurements are 
important for the Standard Model prediction of , which 
differs by more than three standard deviations from a re- 
cent direct BNL measurement [2]. 
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The study of initial-state radiation (ISR) events at an 
e + e~ collider with a fixed center-of-mass (CM) energy 
Ecm allows high-precision measurements of the cross 
section of exclusive hadronic channels for energies below 
the nominal Ecm and is complementary to studies based 
on an energy scan. Use of the ISR technique is discussed 
in Refs. [3-6]. Previously, we used this technique to in- 
vestigate low-multiplicity hadronic processes at effective 
CM energies below 5 GeV [7-15]. 

The ISR cross section cr/. 7 for a specific final state / 
depends on the non-radiative cross section 0/ and is given 
by [5]: 

d<jf n (s,M had ) 2M had 

' = W(s,x,C )-a f (M had ), (1) 

alVlhad S 

where x = 2-E*/y/s, y/s is the nominal CM energy, E* is 

the energy of the ISR photon, 1 and Mhad = ^/s(l — x). 
The radiator function W(s,x, C*), which describes the 



refers to the nominal CM frame. 
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probability for the emission of an ISR photon in the po- 
lar angle range | cos 6*\ < C* , is determined to next-to- 
leading order (NLO) ISR with the PHOKHARA software 
package [16]. 9* is the angle of the emitted 7 with respect 
to the direction of the incoming e~ . Effects of final-state 
radiation (FSR), which are neglected in cq. (1), are stud- 
ied and corrected for as explained in section VI. 

Recently, cross section measurements with approxi- 
mately 1% precision have been presented for the 7r + 7r~ 
final state [15, 17, 18]. These measurements have led 
to a significant reduction in the uncertainty of a^ ad [19]. 
As a consequence, a large relative contribution to the 
uncertainty of aj^ ad now arises from the energy range 
lGeV < E < 2GcV [20]. In this energy range, the 
hadronic cross section is dominated by the exclusive 
e + e~ — > 7T + 7T~7r + 7r~ and 7r + 7r~7r°7r° channels. 

This paper reports results from an ISR analysis of the 
e+e~ — »■ 7r + 7r~7r + 7r~7 process. The four-pion mass 
serves as the effective hadronic Ecu value. It is an up- 
date of our earlier study [8] , based on a data sample that 
is five times larger. We perform a more detailed study 
of systematic effects and thereby obtain significant im- 
provements in both the statistical and systematic uncer- 
tainties. 

The outline of this paper is as follows. Section II de- 
scribes the BABAR detector and the data set used in this 
analysis. The primary event selection is presented in sec- 
tion III. After discussing the background suppression 
(section IV) and the acceptance and efficiency studies 
(section V) , the extraction of the non-radiative cross sec- 
tion is described in section VI. A qualitative analysis of 
the intermediate subsystems and a quantitative measure- 
ment of the J/ip and ip{2S) branching fractions follows 
in section VII. A summary is given in section VIII. 



II. THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with the 
BABAR detector at the PEP-II asymmetric energy e + e~ 
storage rings at the SLAC National Accelerator Labora- 
tory. A total integrated luminosity of 454.3 fb _1 is used, 
comprised of 413.1 fb -1 collected at the T(4S) resonance 
peak, and 41.2 fb~ collected 40MeV below the peak. 

The BABAR detector is described in detail else- 
where [21]. The reconstruction of charged- particle tracks 
is performed with the tracking system, which is com- 
prised of a five-layer silicon vertex tracker (SVT) and 
a 40-layer drift chamber (DCH), both in a 1.5 T axial 
magnetic field. Separation of electrons, protons, charged 
pions, and charged kaons is achieved using Cherenkov 
angles measured with the detector of internally reflected 
Cherenkov light (DIRC) in combination with specific ion- 
ization dE/dx measurements from the SVT and DCH. 
The CsI(Tl) electromagnetic calorimeter (EMC) mea- 
sures the energy of photons and electrons. Muon identi- 
fication is provided by the instrumented flux return. 

A simulation package developed for radiative pro- 



cesses, AFKQED, is used to determine detector accep- 
tance and reconstruction efficiencies. Hadronic final 
states, including 7r + 7r~7r + 7r~7, are simulated based on an 
approach of Czyz and Kiihn [22]. The underlying model 
assumes dominance of the ai(1260)7r final state as was 
reported in Refs. [23, 24]. Due to the dominant decay 
ai(1260) — > p°n, each event contains one pair of pions 
from p° decay. The simulation of multiple soft-photon 
emission from the initial state is performed via a struc- 
ture function technique [25, 26]. Extra radiation from 
the final state particles is simulated by the PHOTOS 
package [27]. The accuracy of the radiative corrections 
is about 1%. A new version of PHOKHARA [16], which 
incorporates the results of recent studies [8] on interme- 
diate resonances, is used to investigate the influence of 
these intermediate resonances on the acceptance. 

The simulated events are subjected to simulation of 
the detector [28] and the same analysis procedures as the 
data. Variations in detector and background conditions 
over the course of the experiment are modeled. 

A large number of potential background processes are 
simulated, including the ISR processes K + K ~n + -K~j 
and K®K ± 7r =F 7. Other ISR background channels are also 
examined. Either the remaining number of events after 
the event selection is negligible in comparison to other 
uncertainties, or dedicated methods for background re- 
jection arc implemented as described in section IV. Non- 
ISR backgrounds resulting from e+e~ — >• qq (q =u,d,s,c) 
continuum events are modeled using the JETSET gener- 
ator [29], while those from e+e~ — > t+t~ arc modeled 
with KORALB [30]. The cross sections of the ISR chan- 
nels and the branching fractions of the inclusive processes 
leading to final states similar to our signal are known 
with about 10% accuracy or better, which is sufficient 
for the purpose of this measurement. The contribution 
from T(45) — > 7r + 7r~7r + 7r~7 events is negligible. 



III. PRIMARY EVENT SELECTION AND 
KINEMATIC FIT 

Charged tracks are selected by requiring that they orig- 
inate from the collision region (transverse distance of 
closest approach to the nominal interaction point dr < 
1.5 cm, and in the beam direction dz < 2.5 cm), and that 
they have a polar angle 9 c h in the well-understood ac- 
ceptance region of the detector (0.5rad < 6* c h < 2.4rad). 
The coordinate system has the z-axis in the direction 
of the incoming e~ beam. Tracks with transverse mo- 
menta less than lOOMeV/c or that are consistent with 
being an electron are rejected. Photon candidates are 
required to have a minimum energy E 1 ^cm > 50MeV. 
The ISR photon candidate is restricted to the polar an- 
gle range inside the well-understood acceptance region 
of the EMC (0.35 rad < 9 1 < 2.4 rad) and a minimum 
energy of £?isr > 3 GeV is required. 

Radiative Bhabha events are suppressed by requiring 
the two most energetic tracks of the event not to be 
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identified as electrons. The minimum angle between a 
charged track and the ISR photon Atp is required to sat- 
isfy Atp > 1.0 rad, in order to select the back-to-back 
topology between the charged tracks and photon typical 
for ISR events. 

A kinematic fit procedure with four constraints (4C) is 
applied to events with four tracks satisfying these crite- 
ria. The constrained fit uses the measured momenta and 
directions of the charged particles and the ISR photon 
and the corresponding error matrix to solve the energy- 
momentum equation. To provide accurate photon pa- 
rameters to the kinematic fit, a precise alignment and 
an energy calibration of the EMC are performed using a 
/i + /i~7 sample. This improves the data-simulation agree- 
ment in the goodness-of-fit distributions (see below) of 
the kinematic fit. The /i + /i~7 sample is also used to 
identify and measure differences between the data and 
the simulation (MC) in the photon detection efficiency. 

The kinematic fit is performed assuming the 
7T + 7r~7r + 7r~7 signal hypothesis. If two tracks are iden- 
tified as kaons the fit is also performed under the 
K + K~tt + it~"/ hypothesis. The fitting routine returns 
a goodness-of-fit quantity xin an d X2K2t^ respectively, 
which we use to select signal events and to suppress 
K + K~Tr + Tr~ r y background. 

Figures 1(a) and (b) show the measured xl^ distri- 
bution before background subtraction in comparison to 
the distribution from the simulated signal sample. The 
results are shown on logarithmic and linear scales. The 
distributions are normalized to the number of events with 
X4 ff < 10. Figure 1(c) shows the x\-k distribution for the 
main backgrounds. The difference between measured and 
simulated samples is shown in Fig. 1(d) for > 10- F° r 
large values of xl^, the difference is approximately flat, 
and is consistent with the sum of backgrounds shown in 
Fig. 1(c). At small xin values, resolution effects cause 
the difference to decrease. This is studied with a clean 
sample of four pions and one photon where the resolution 
effect is visible for x%k < 20. To avoid a bias due to the 
resolution difference between the data and simulation, 
the requirement x%n < 30 is used in the event analysis. 
In addition, background channels are suppressed using 
dedicated vetoes and particle identification (PID) selec- 
tors, as discussed in section IV. 



IV. BACKGROUND REJECTION 

Due to large variations in the signal-to-background ra- 
tio across different Air mass regions, different strate- 
gies are implemented to eliminate background, depend- 
ing on the region. Figure 2(a) shows the dis- 
tribution under the 7r + 7r~7r + 7r~7 hypothesis with the re- 
quirement x\ix < 30. The results are shown for the data 
and for the sum of the simulated ISR channels e + e~ — » 
K + K~tt + it~j, e + e~ — > K^K^tt^j, e + e~ — > ujn + 7r~j 
and the non-ISR continuum background. In Fig. 2(b) 
these background contributions are shown individually. 




FIG. 1. (a,b) The x\tt distribution for data without back- 
ground subtraction (triangles) and signal MC (histogram) on 
a logarithmic and linear scale, (c) MC distributions of the 
principal backgrounds, (d) The difference between the mea- 
sured and the simulated distributions from part (a,b). The 
dashed and solid vertical lines indicate the boundaries with 
the < xln < 10 region used for normalization and the 
x\tt < 30 requirement, respectively. 



The xl-n- requirement effectively eliminates the W7r + 7r~7 
background: according to simulation, only 450 W7r + 7r~7 
events remain after applying this restriction, in compar- 
ison to 260 000 signal events with a very similar 
distribution, leading to a relative contribution of less 
than « 0.2%. At very low invariant masses, < 
l.lGeV/c 2 , the background contributions from 7r + 7r~77 
and 7r + 7r~7r°7 are large and not included in Fig. 2. In 
this threshold region, the signal-to-background ratio is 
approximately 1:5. In the peak region, l.lGeV/c 2 < 
Af^ < 2.2GeV/c 2 , the background contamination is at 
the level of 3-4%, dominated by e + e~ — > K + K~n + it~^ 
and e + e~ — » K^K^n^j. At high M^, the background 
level rises to about 10%, mainly due to the additional 
contribution of uds-continuum events. 

A clear peak from J/tp — > ir + TT~ir + n~ is visible in the 
distributions of Fig 2. In addition, the data exhibit a 
narrow peak in the ip(2S) mass region, due to the decay 
ip(2S) — > 7r + 7r~ J/tp with J/ip — > /i + and misidcntificd 
muons. The rejection of this background is described in 
subsection IV B. 
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FIG. 2. (a) The tt + -k~tt + -k~ invariant mass distribution of 
7T + 7r _ 7r + 7r^7 ISR events (triangles) for the main selection and 
the sum of simulated K + K~tt + tt~ 7, KgK if^y, oj7r + 7r _ 7 
and non-radiative uds-continuum background (histogram). 
The large 7r + 7r~77 and 7r + 7r~7r°7 background contribution 
in the threshold region (M^ < l.lGeV/c 2 ) is not included 
in the simulation, (b) The individual simulated background 
channel contributions from top to bottom as indicated in the 
legend. 



A. Background in the peak region 

(l.lGeV/c 2 < M 4T < 2.2GeV/c 2 ) 

The dominant background in the peak region is from 
the ISR processes e + e~ — > K + K~Tr + ir~~f and e + e~ — > 
K®K ± ir T '-f. We utilize two different approaches to eval- 
uate the background and use the difference between the 
methods to estimate the systematic uncertainty related 
to the background subtraction. The background subtrac- 
tion itself is performed using method 2. 



Method 1: subtract simulated background 

For background subtraction via method 1, we use re- 
cent measurements of the most important background 
channels, K+K~ tt+tt^ [13] and K°K ± ir ¥ 'y [14]. We 
tune AFKQED according to these measurements and use 
the resulting predictions to evaluate the contributions of 
these channels to the spectrum. 



Method 2: background suppression 

This method is a hybrid approach. We impose specific 
requirements in order to suppress the A" + A" _ 7r + 7r _ 7 and 
K®K ir^j backgrounds: the so-called A" + A"~7r + 7r~7 
and A"°A" ± 7r =F 7 vetoes. The remaining background is 
then subtracted according to method 1. 

A"gA" ± 7r =F 7 events are vetoed if one of the tracks is 
identified as a K^. The identification algorithm has an 
efficiency of 80-90% per track and a it mis-identification 
probability of 0.25%. In addition we require the invariant 
mass formed from two of the three remaining tracks to 



lie within 35 MeV/c 2 of the nominal mass [31]. Ac- 
cording to simulation, this method removes 74% of the 
K^K ± tt t j background but only 1% of signal. 

K + K ~7r + 7r~7 events are rejected if two oppositely- 
charged tracks arc identified as kaons and the require- 
ment X2K2tt < 10 is fulfilled. The identification algo- 
rithm has an identification efficiency of 85-95% per track 
and a 7r mis-identification probability of 1%. This re- 
quirement removes less than 0.1% of the signal, but 55% 
of the A" + A"~7r + 7r~7 background. 

Comparing the results of the two methods yields a sys- 
tematic uncertainty of 1.0% on the e + e~ — > 7r + 7r~7r + 7r~ 
cross section in the peak region. 



B. Background at large invariant masses 

(M i7! > 2.2GeV/c 2 ) 

From Fig. 2(b) it is seen that the uds continuum 
background is significant in the > 2.2GeV/c 2 

mass region. The largest contribution is from e + e~ — > 
7r + 7r~7T + 7r~7r events in which one of the photons from 
the 7T° — > 77 decay is mistaken for the Jisr- The simi- 
lar kinematic configuration of the uds-continuum events 
causes a peak at small values of x%k- We estimate this 
background by measuring the ir° yield from a fit to the 
7T° mass peak in the two-photon invariant mass distribu- 
tion of the ISR photon candidate with any other photon 
candidate. We then scale the MC 7r° rate to give the 
same yield as observed in the data. The corresponding 
scaling factor is extracted with a relative uncertainty of 
about 22%. The normalized wds-continuum MC sam- 
ple is used for background subtraction. The correspond- 
ing systematic uncertainty on the cross section measure- 
ment is 0.5% in the peak region and rises up to 1.5% for 
M 47r > 2.8GeV/c 2 . 

There is an additional peaking background contribu- 
tion in the charmonium region. Figure 3(a) shows the in- 
variant mass distribution under the 7r + 7r~7r + 7r~ hypoth- 
esis for the data. Peaks corresponding to the J/0 and 
0(25) resonances are clearly seen. In Fig. 3(b) the invari- 
ant 7r + 7r~ mass distribution for events within the 
0(25) (3.65GeV/c 2 < M 47r < 3.75GeV/c 2 ) mass region 
are plotted. The distribution contains four entries per 
event. The J/0 peak is clearly visible. The J/0 — > tc + it~ 
branching fraction is about a factor 400 smaller than that 
of J/0 — > [31]. No PID selector is used in the 

sample selection. In addition, the kinematic fit does not 
suppress 7r + 7r~^ + /i _ events, because the mass difference 
between muons and pions is only 34MeV/c 2 , which is 
negligible due to the fact that they are highly relativis- 
tic, being emitted from the J/0 . Thus, the observed peak 
is likely dominated by 0(25) — > n + ir~ /i~ decays. To 
remove this background, we reject events that have 
in the 0(25) region and an value consistent with 
the J/0 . The corresponding distribution after the 
0(25) veto is shown in Fig. 3(c). 
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FIG. 3. (a) Invariant mass distribution under the 7r + 7r~7r + 7r~ 
hypothesis for the data; (b) tv + it~ invariant mass distribution 
for events with in the tp{2S) mass region represented by 
the vertical lines in (a) (4 entries per event); (c) Same as part 
(a) after removing entries of the invariant 7r + 7r _ -mass in the 
J/%j) mass region represented by the vertical lines in (b). 



the track. The first requirement is not sufficient to re- 
ject background events with highly-energetic electrons or 
positrons, which have a vertex probability similar to non- 
conversion events. MC study shows that the combination 
of both requirements yields a background rejection larger 
than 99% while removing less than 6% of signal. 



Method 2: Pion identification 

The second method requires all four tracks to be 
identified as pions, using a selector with an efficiency 
of 97-99% and an electron mis-identification probability 
of 5-7% depending on the transverse momentum of 
the track. There is, however, a difference in efficiency 
between the data and MC simulation of approximately 
0.5 - I % per track (2 - 4% shift per event). Therefore 
this selector is not used in the peak region, where the 
background contribution is very low. 

Both methods remove a large fraction of the conver- 
sions and yield results that are consistent with each other. 
We present our primary results using the pion selection 
method, as the systematic uncertainties on its inputs are 
better understood. Comparing the cross section results 
using the two methods leads to a systematic uncertainty 
estimate of 3%. It should be noted that the p(770)° peak 
is strongly suppressed, but is still visible as a shoulder 
after background subtraction, as shown in Fig. 4. 



C. Background at small invariant masses 

(M 4n < l.lGeV/c 2 ) 

Background in the threshold region, i.e., below 
I.I GeV/c 2 , is dominated by two processes: 7 conversion 
of a real photon in the detector material, and conver- 
sion of a virtual photon at the primary interaction ver- 
tex (Dalitz conversion). In e + e~ — > 7r + 7r~7r°7isR events, 
one of the photons from the ir° decay can convert in 
the detector material into an e + e~ pair, and both the 
e + and the e~ can be misidentified as pions. Moreover 
e + e _ — > 7r + 7r~7_r,s_R7 events in which the non-ISR pho- 
ton converts into an e + e~ pair contribute a similar back- 
ground. The second source of background in this mass 
region is e + e~ — > ir + Tt~"fisRe + e~ events in which the 
e + e~ pair arises from a Dalitz conversion process. Two 
methods are used to remove these background channels. 



Method 1: Vetoes for Dalitz conversion and conversion in 
the detector 



The first method vetoes events with either a primary 
vertex probability less than 10~ 8 or with two identified 
electrons. Electrons are identified using an algorithm 
with an efficiency of 99% and a pion mis-identification 
rate of 5-10% depending on the transverse momentum of 
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FIG. 4. 7r + 7r~7r + 7r~7 cross section in the threshold region 
for the data (squares) and for the data after requiring each 
track to satisfy a dedicated pion selection algorithm (trian- 
gles), corresponding to the method 2 selection criterion. The 
indicated uncertainties are statistical. 
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FIG. 5. (a) The \ 2 distribution for data after background 
subtraction (points), signal MC (histogram) and the sum of 
simulated background channels (hatched histogram), (b) The 
difference between data with background subtraction and sig- 
nal MC. The dashed line indicates the boundary with the 
< xlir < 10 region used for normalization. 
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D. Background subtraction summary 

Figure 5(a) shows the xl-n distribution both for the 
data after background subtraction and the signal simula- 
tion. Compared to Fig. 1(a), the difference between the 
data and the signal simulation is reduced. There is still 
some remaining background at larger xl-x values as can 
be seen in Fig. 5(b). This is consistent with the uncer- 
tainties in the background-subtraction methods as previ- 
ously described. To be conservative, we assume that the 
remaining background is uniform as a function of xItt I n 
this case, we find that the associated uncertainty on the 
cross section is less than 0.4% for < 2.8 GeV/c 2 and 
4.0% for 2.8GeV/c 2 < M 47I < 4.5GeV/c 2 . 



V. ACCEPTANCE AND EFFICIENCIES 



FIG. 6. The efficiency as a function of the n + ir~ir ir~ invari- 
ant mass. The vertical line at — l.lGeV/c 2 separates 
the regions with and without a pion PID requirement for the 
reconstructed tracks. The curves on either side of the vertical 
line are fits based on the sum of a Gaussian and a polynomial. 

As mentioned in section II, it is assumed that the 
7r + 7r~7r + 7r~7 hadronic final state arises from the decay 
of various intermediate resonances such as ai(1260)7r. 
The relative contributions of intermediate states are dis- 
cussed in section VII. Different intermediate states 
might exhibit different results for the angular distribu- 
tions of final-state particles. The limited knowledge of 
the hadronic substructure hence corresponds to a sys- 
tematic uncertainty in the evaluation of the global ef- 
ficiency. This effect is estimated through comparison 
of the predictions from AFKQED and a new version of 
PHOKHARA [16], which contain different intermediate 
resonances. The effect is observed to be negligible. 



The full chain of experimental requirements is applied 
to the signal MC sample. We define the ratio of the num- 
ber of selected events divided by the number of events 
without applying any requirements as the global effi- 
ciency. Figure 6 shows the global efficiency as a function 
of determined with the simulation. The decrease of 
efficiency in the threshold region, which corresponds to 
the highest energy ISR photons, is due to the fact that 
the four tracks of the hadronic system recoil in a narrow 
cone opposite to the direction of the ISR photon. Since 
ISR photons are preferentially emitted at small polar an- 
gles, the efficiency for detecting all four tracks within 
the fiducial volume of the detector decreases accordingly. 
The decrease in the global efficiency at large values of 
M^-x can be explained with similar arguments, because 
the opening angle of the hadronic system increases with 
decreasing ISR photon energy, increasing the probabil- 
ity of losing one of the four tracks. The discontinuity at 
Matt = 1.1 GeV/ c 2 is due to the pion PID requirement for 



A. Photon efficiency 

A dedicated study is performed to determine the pho- 
ton efficiency in data and simulation. Detector ineffi- 
ciencies due to inactive material between crystals, non- 
functioning crystals, and conversions in inner detector 
structures are investigated. This study is performed with 
a /i + /J,~"fisn sample, in which an identified ISR photon is 
not required in the event selection. A kinematic fit with 
one constraint (1C) based on the kinematic information 
of the charged tracks is performed in order to calculate 
the energy and direction of the photon. This fit predic- 
tion is compared to the measured photon information to 
extract the inefficiency of the photon reconstruction. In 
Fig. 7 the photon inefficiency is shown as a function of 
the polar angle of the photon for the data and simula- 
tion. The inefficiency curve is not smooth due to the 
effect of the inactive regions between EMC crystal rings. 
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Some gaps between rings, which are visible in the data 
as peaks of high inefficiency, are not properly simulated, 
especially in the forward region of the detector. As a 
function of M 47r , we observe a uniform inefficiency differ- 
ence between the data and the simulation with an average 
value of Ae 7 = (1.34 ± 0.03 sfQt ± 0.37 syst)%- 
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FIG. 7. Photon inefficiency as a function of the ISR photon 
polar angle for the data (circles) and signal MC (crosses). The 
binsize is increased to 6mrad for ^photon > 2.1 rad. 



B. Tracking efficiency 

Track reconstruction and effects like nuclear interac- 
tions and energy loss of tracks traversing the detector 
volume are not simulated perfectly. As a consequence, 
the tracking efficiency is slightly different for the data 
and MC. This difference is investigated with a dedicated 



study of ISR e H 



events with one miss- 



ing track. The missing momentum and direction of the 
lost pion is calculated using a constrained kinematic fit 
with one remaining constraint (1C). The angular and mo- 
mentum dependent distribution of the inefficient events 
is extracted. 

Small differences in tracking efficiency between the 
data and simulation arc due to an imperfect description 
of track loss when tracks overlap in azimuth. The aver- 
age difference between the data and MC is uniform as 
a function of the transverse momentum and the polar 
angle of the tracks and is determined on average to be 
Ae trfc = (0.75 ± 0.05 stQt ± 0.34 syst )%. 



VI. CROSS SECTION 

As described in section I, the non-radiative cross sec- 
tion is related to the radiative cross section according 
to: 

s 1 



2Ma 



W{s,x,C*) 



(2) 



The radiative cross section is: 
dcr 4 ^, 7 (M 47r ) _ dN^ n (M 4n ) 



1 



dM. 



4 71 



dM 4n 



(1 



Srad,FSR) 

where dN^.-f is the number of selected events, e the 
global efficiency corrected for tracking and photon ef- 
ficiency differences between the data and MC, and 
Srad,FSR the radiative corrections including LO- and 
NLO-FSR effects. The MC event generator interfaced 
with the detector simulation is based on a modified ver- 
sion of the EVA code [22]. It contains collinear NLO-ISR 
corrections based on the structure function technique, as 
well as FSR corrections based on PHOTOS [27]. Radia- 
tive corrections due to NLO-ISR on the radiator func- 
tion obtained with this event generator are compared to 
the PHOKHARA generator [16], which includes the full 
NLO-ISR corrections. With our selection, a (1.0 ±0.2)% 
difference of the radiator function between the two gen- 
erators is observed. We apply a correction to account for 
this difference. 
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FIG. 8. Ratio of the simulated radiative cross section with 
only the ISR contribution to the cross section with both ISR 
and FSR. 

Figure 8 shows the ratio of the simulated ISR radiative 
cross section and the cross section including additional 
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FSR (PHOTOS). The FSR leads to a shift of events to- 
wards lower invariant masses in the radiative cross sec- 
tion, due to the fact that the measured invariant mass is 
smaller than the effective Ecm for events with FSR. 

After applying all radiative corrections and account- 
ing for the relevant differences in efficiencies between the 
data and simulation, we apply eq. (2) and extract the 
non-radiative e + e~ — > 7r + 7r~7r + 7r~ cross section. The 
result is shown in Fig. 9. The measured cross section 
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cross section. The uncertainties associated with back- 
ground subtraction are discussed in section IV. The 3.0% 
tracking efficiency difference between data and MC has 
an uncertainty of 1.4%. The photon efficiency correction 
is 1.3 ± 0.4%. The total luminosity is measured with a 
precision of 1.0%. A 1.0 ±0.2% difference is observed be- 
tween the radiator functions computed with AFKQED 
and PHOKHARA. The effect of additional FSR is esti- 
mated using PHOTOS, resulting in the correction shown 
in Fig. 8 and a systematic uncertainty of 0.5%. The re- 
quirement xln < 30 leads to a systematic uncertainty 
of 0.3%. The uncertainty on the global efficiency is esti- 
mated to be 1.0% in the central region, increasing to 10% 
in the low mass region < 1.1 GeV/c 2 due to an ob- 
served efficiency decrease of up to 10%. A conservative 
uncertainty of 10% to account for the total acceptance 
decrease in this region is also applied. 

Assuming no correlation between the various contribu- 
tions to the systematic uncertainty of the cross section, 
its total is found to be 10.7% for < l.lGeV/c 2 , 

2.4% for l.lGeV/c 2 < M 4w < 2.8GeV/c 2 , 5.5% for 
2.8GeV/c 2 < M 47T < 4.0GeV/c 2 and 8.5% for higher 
invariant masses. Individual contributions to the sys- 
tematic uncertainties contribute in a correlated way on 
the whole mass range, with the exception of the global 
analysis efficiency, for which it does not. Therefore for 
M 47r > l.lGeV/c 2 a 100% correlation can be assumed, 
while for M 4v < l.lGeV/c 2 where the global efficiency 
dominates, it can be assumed to be uncorrelated. 

B. Comparison with the existing e + e~ data 



FIG. 9. The Ecm dependence of the dressed e + e~ — ¥ 
tt + tv~tt + -k~ cross section measured from the ISR data. The 
uncertainties are statistical. 

includes the contributions of vacuum polarization. The 
cross section used in the dispersion integral for a M does 
not include vacuum polarization and thus we need to 
apply a correction. We define the undressed cross sec- 
tion & 4 ^(Ecm) by correcting the measured or dressed 
cross section ^^{Ecm) for vacuum polarization effects 
according to: 

ud (F \ ( a (°) V_ <r(ffgM) 

OaAEcm) - u^yhjcM) ■ — tt; r - t tt; r 

\a(hcM)J o vac {EcM) 

(4) 

where a(EcAt) is the electroweak coupling strength 
at Ecm- The correction due to vacuum polarization 
dvac(EcM), which can be found in Ref. [32], is applied. 
Our results for the dressed and the undressed cross sec- 
tions are presented in Table I. 

A. Systematic corrections and uncertainties 

Table II presents the complete list of corrections and 
systematic uncertainties that are included in the dressed 



In Fig. 10 the extracted non-radiative <r(e + e~ — > 
7r + 7r~7r + 7r~) cross section is shown, in comparison with 
the previous BABAR result [8] and the results from fixed- 
energy e + e _ experiments. Our results agree within the 
uncertainties with our previous measurement, which they 
supersede. Our results are consistent with and higher 
in precision than the direct e + e~ cross section mea- 
surements made at VEPP-2M by OLYA [33], ND [34], 
SND [35], CMD [36], and CMD-2 [37-39], at DCI by 
M3N [40], DM1 [41], and DM2 [42] and at Adone by 
GG2 [43]. 

C. Influence on the prediction of a M 

Using the result for the e + e~ — > 7r + 7r~7r + 7r~ cross 
section obtained in the present study, we compute the 
contribution of this channel to the anomalous magnetic 
moment of the muon via a dispersion relation us- 
ing the HVPTool program [44] in the energy region 
0.6 GeV < E CM < 1.8 GeV. We find: 

ajTV+TT-TT+O = (13.64 ±0.03 stat ±0.36 syst ) x lO" 10 . 

Our result is more precise than the current world average 
for this quantity: (13.35±0.10 stat ±0.52 S j /st ) x 10~ 10 [19], 
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TABLE I. Summary of e + e —¥ tt + it tt + it cross section measurement. Dressed (with VP) and undressed (without VP) cross 
sections are presented with statistical uncertainties only. 
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1.95 ± 0.07 


3812.5 


0.19 ± 0.02 


0.18 ± 0.02 


1237.5 


12.62 ± 0.20 


12.17 ± 0.20 


2537.5 


1.87 ±0.07 


1.80 ± 0.06 


3837.5 


0.18 ± 0.02 


0.17 ± 0.02 


1262.5 


14.56 ± 0.22 


14.05 ± 0.21 


2562.5 


1.71 ± 0.06 


1.65 ± 0.06 


3862.5 


0.18 ± 0.02 


0.17 ± 0.02 


1287.5 


16.39 ± 0.23 


15.83 ± 0.22 


2587.5 


1.85 ± 0.06 


1.77 ± 0.06 


3887.5 


0.21 ± 0.02 


0.20 ± 0.02 


1312.5 


19.06 ± 0.25 


18.41 ± 0.24 


2612.5 


1.79 ± 0.06 


1.72 ± 0.06 


3912.5 


0.20 ± 0.02 


0.19 ± 0.02 


1337.5 


21.14 ± 0.26 


20.42 ± 0.25 


2637.5 


1.62 ±0.06 


1.56 ± 0.06 


3937.5 


0.15 ±0.02 


0.14 ±0.02 


1362.5 


23.37 ± 0.27 


22.59 ± 0.26 


2662.5 


1.43 ± 0.06 


1.37 ± 0.05 


3962.5 


0.14 ± 0.02 


0.14 ± 0.01 


1387.5 


25.76 ±0.28 


24.90 ± 0.28 


2687.5 


1.31 ±0.05 


1.26 ± 0.05 


3987.5 


0.16 ±0.02 


0.16 ± 0.02 


1412.5 


27.53 ± 0.29 


26.61 ± 0.29 


2712.5 


1.30 ± 0.05 


1.26 ± 0.05 


4012.5 


0.17 ± 0.02 


0.16 ± 0.02 


1437.5 


29.95 ± 0.30 


28.96 ± 0.30 


2737.5 


1.21 ± 0.05 


1.16 ± 0.05 


4037.5 


0.12 ± 0.01 


0.11 ± 0.01 


1462.5 


30.32 ± 0.31 


29.32 ± 0.30 


2762.5 


1.17 ±0.05 


1.13 ± 0.05 


4062.5 


0.20 ± 0.02 


0.19 ± 0.02 


1487.5 


32.04 ± 0.31 


30.97 ± 0.30 


2787.5 


1.17 ±0.05 


1.12 ± 0.05 


4087.5 


0.13 ± 0.01 


0.12 ± 0.01 


1512.5 


30.98 ± 0.31 


29.93 ± 0.30 


2812.5 


1.09 ±0.05 


1.05 ± 0.05 


4112.5 


0.14 ±0.02 


0.13 ± 0.01 


1537.5 


30.11 ± 0.30 


29.06 ± 0.29 


2837.5 


1.07 ±0.05 


1.04 ±0.05 


4137.5 


0.14 ± 0.02 


0.13 ± 0.01 


1562.5 


28.26 ± 0.29 


27.26 ± 0.28 


2862.5 


0.96 ± 0.05 


0.93 ± 0.04 


4162.5 


0.14 ± 0.01 


0.13 ± 0.01 


1587.5 


26.81 ± 0.28 


25.86 ± 0.27 


2887.5 


0.89 ±0.04 


0.86 ± 0.04 


4187.5 


0.15 ±0.02 


0.14 ± 0.01 
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1637.5 


22.69 ±0.26 


21.89 ±0.25 


2937.5 


0.88 ±0.04 


0.85 ± 0.04 


4237.5 


0.13 ±0.01 


0.12 ±0.01 


1662.5 


20.95 ± 0.25 


20.19 ± 0.24 


2962.5 


0.77 ±0.04 


0.75 ± 0.04 


4262.5 


0.13 ±0.01 


0.12 ±0.01 


1687.5 


18.78 ±0.23 


18.09 ± 0.22 


2987.5 


0.82 ±0.04 


0.81 ±0.04 


4287.5 


0.13 ±0.01 


0.12 ±0.01 


1712.5 


17.25 ±0.22 


16.61 ± 0.21 


3012.5 


0.75 ±0.04 


0.74 ±0.04 


4312.5 


0.11 ±0.01 


0.11 ± 0.01 


1737.5 


15.33 ±0.21 


14.75 ± 0.20 


3037.5 


0.71 ±0.04 


0.71 ±0.04 


4337.5 


0.11 ±0.01 


0.11 ± 0.01 


1762.5 


13.37 ±0.19 


12.86 ± 0.19 


3062.5 


0.62 ±0.04 


0.66 ± 0.04 


4362.5 


0.09 ±0.01 


0.09 ± 0.01 


1787.5 


11.61 ±0.18 


11.17 ±0.17 


3087.5 


1.93 ±0.06 


2.30 ± 0.08 


4387.5 


0.10 ±0.01 


0.10 ±0.01 


1812.5 


10.23 ±0.17 


9.84 ±0.16 


3112.5 


1.30 ±0.05 


1.03 ±0.04 


4412.5 


0.11 ±0.01 


0.10 ±0.01 


1837.5 


8.87 ±0.15 


8.53 ±0.15 


3137.5 


0.62 ± 0.04 


0.55 ±0.03 


4437.5 


0.09 ±0.01 


0.09 ± 0.01 


1862.5 


7.67 ±0.14 


7.37 ±0.14 


3162.5 


0.59 ± 0.03 


0.54 ±0.03 


4462.5 


0.09 ±0.01 


0.08 ± 0.01 


1887.5 


7.29 ±0.14 


7.02 ±0.13 


3187.5 


0.51 ±0.03 


0.47 ±0.03 


4487.5 


0.10 ±0.01 


0.09 ± 0.01 
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TABLE II. Summary of systematic corrections and uncertainties in per centage. 



M 47I < l.lGeV/c 2 < 2.8GeV/c 2 < 4.0GeV/e 2 < 4.5GeV/e 2 



K + K~TV + 7V'J, _K"°A" ± 7T T 7 


± 1.0 


± 1.0 


± 3.0 


± 7.0 


continuum bkg 




± 0.5 


± 1.0 


± 1.5 




± 3.0 








additional bkg 


± 0.4 


± 0.4 


± 4.0 


± 4.0 


tracking efficiency 


+3.0 ± 1.4 


+3.0 ± 1.4 


+3.0 ± 1.4 


+3.0 ± 1.4 


photon efficiency 


+1.3 ± 0.4 


+1.3 ± 0.4 


+1.3 ± 0.4 


+1.3 ± 0.4 


£ 


± 1.0 


± 1.0 


± 1.0 


± 1.0 


AFK-PHOK-difference 


-1.0 ± 0.2 


-1.0 ± 0.2 


-1.0 ± 0.2 


-1.0 ± 0.2 


FSR corrections 


± 0.5 


± 0.2 


± 0.1 


± 0.1 


Xlr < 30 


± 0.3 


± 0.3 


± 0.3 


± 0.3 


global efficiency 


± 10.0 


± 1.0 


± 1.0 


± 1.0 


sum 


± 10.7 


± 2.4 


± 5.5 


± 8.5 
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FIG. 10. The e e~ — > it ty~tt tt~ cross section in compar- 
ison to other experiments. Only statistical uncertainties are 
shown. 
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FIG. 11. Invariant tt + it~tt ± and ty + tt~ mass combinations vs. 
invariant 7r + 7r~7r + 7r - mass for the data without background 
subtraction (left) and signal MC (right). 



where the first uncertainty is statistical and the second 
systematic. 



VII. INVARIANT MASSES AND 
CHARMONIUM BRANCHING RATIOS 

Different invariant mass combinations have been stud- 
ied in the data and MC simulation to search for states 
not included in the MC model. In the following, we 
present a general qualitative search for these hadronic 
structures. We then consider a more detailed study of 



the J/ip and ip(2S) background subtraction and efficiency 
corrections. Finally, we determine the branching frac- 
tions Bjfy^+K-w+w- and Bi>(2S)^J/ipn+n- and perform 
a scan for additional resonances at high invariant masses. 



A. Substructures 

The scatter plots in Fig. 11 display distributions of the 
invariant 7r + 7r~7r ± and tt + tt~ masses versus the invariant 
ir + ir~ir + iT~ mass for the data and MC. The p(770)° band 
is clearly visible in the tt + tt~ mass distribution of the 
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data and MC. In general, good agreement is seen except 
for the J/tp decay, which is not simulated. 

In a more detailed study the 7r + 7r~7r + 7r~ mass spec- 
trum is divided into five intervals: 

1. 1.0-1.4 GcV/c 2 : low mass region 

2. 1.4-1.8 GeV/c 2 : peak region of the cross section 

3. 1.8-2.3 GeV/c 2 : high mass shoulder 

4. 2.3-3.0 GeV/c 2 

5. 3.0-4.5 GeV/c 2 : without the narrow region around 

w 

Figure 12 shows the one-dimensional distributions 
from the five regions for the two- and three-pion invariant 
masses in comparison with MC [22]. 

The 7T + 7T~n ± invariant mass distribution is shown in 
the leftmost column of Fig. 12. In the low mass region, 
1.0GeV/c 2 < M i7r < 1.4GcV/c 2 , there is not enough en- 
ergy to allow production of the ai(1260) ± . At higher 
M 47r , the contribution of the ai(1260) ± becomes visi- 
ble. It is observed as a peaking structure with mass 
and width M 37r w 1300MeV/c 2 and T w 200 MeV. In 
comparison, the average mass value in the PDG [31] is 
1.230 ± 0.040 GeV/c 2 , with results from individual exper- 
iments that vary between 1.04 and 1.33 GcV/c 2 . The cor- 
responding width varies between 250 and 600 MeV. In 
our simulation, the parameters M = 1.33GeV/c 2 and 
T = 570 MeV are used, which are determined from a 
combined analysis of the CLEO and CMD-2 data [24]. 
Our results seem to favor a lower ai(1260) ± mass and a 
smaller width. 

In the 7r+7r~ invariant mass distributions shown in 
the middle column of Fig. 12, four entries are present 
per event. At low A-n mass and in the peak region 
only a single resonance, the p(770)°, is observed. At 
larger An mass, a second peaking structure appears at 
M 27T ~ 1270 MeV/ c 2 , which most likely corresponds to 
the /2(1270). This resonance is not simulated by our 
MC. It is observed that over the entire An mass range, 
approximately 25% of the entries are in the p(770)° peak. 
p(770)°p(770)° production is not allowed due to C-parity 
conservation, leading to the conclusion that in each event 
one p(770)° meson is present. 

To investigate the possible presence of the 
/ 2 (1270)p(770)° final state, the tt+tt^ combination 
is plotted for the case that there is another 7r + 7r~ 
combination within ±25 MeV/ c 2 of the p(770)° mass, 
745MeV/c 2 < M 27r < 795MeV/c 2 . The results are 
shown in the rightmost column of Fig. 12. An artificial 
dip at M 27r w 770 MeV/ c 2 due to the selection of the 
p(770)° is observed. The /2(1270) resonance is visible 
as a shoulder in the 1.8GeV/c 2 < M 47r < 2.3GeV/c 2 
mass region. It is even more prominent in the 
2.3GeV/c 2 < M 47r < 4.5GeV/c 2 region, where 
the energy is large enough to allow direct pro- 
duction of / 2 (1270)p(770)°. A sharp falloff in the 



M2tt spectrum just below lGeV/c 2 is visible in the 
1.8GeV/c 2 < M 47r < 2.3GeV/c 2 region. This might be 
due to interference with the /o(980) final state. A partial 
wave analysis, which is beyond the scope of this paper, 
would be necessary to determine the structure of the 
individual intermediate states. A qualitative comparison 
with the MC model of Ref. [16] indicates a somewhat 
better agreement with the data than that shown in the 
rightmost column of Fig. 12, apart from an overestimate 
of the contribution of the / (1300)p(770)° final state. 

B. J/tp and ip(2S) 

Figure 13(a) displays the cr(e + e~ — > 7r + 7r~7r + 7r~) cross 
section as a function of the 4tt mass, in the vicinity of the 
J/tp meson. The measured width of the J/tp (w 15 MeV) 
is dominated by the track momentum resolution, the in- 
trinsic width being Tj/^p = 93keV [31]. The small tail 
towards higher masses is mostly from extra radiation, 
which is assigned to the hadronic system by the fit. We 
describe the J/tp peak and the non-resonant 7r+ tt~ ~ir + ~ir~ 
contribution with the sum of two Gaussians and a lin- 
ear term, respectively. This allows the integrated partial 
cross section cr^f = f °° dM^a 3 ^ (M^) and the elec- 
tronic width of J/tp — )• 7r + 7r~7r + 7r~ to be extracted: 

n J/4> 

N(J/tp -> 7r + 7T~7r + 7r~) 
dC/dE ■ e M c 
= (48.9 ± 2.1 stat ± 1.0 syst ) McV/c 2 nb (6) 

and 

^ J/tp — >7T+ 7T~ 7T+7r~ ' ^ ee 

N(J/tf) -> Tr+TT-TT+Tr-) • M 2 ^ 
_ 6tt 2 • dC/dE -6MC-C 
= (20.4±0.9 stat ±0.4 SJ/st )eV. (7) 

For the above, the value Mj^ = 3096.92 ± 
0.01 MeV/c 2 [31] and the conversion constant C — 
3.8938 x 10 11 MeV 2 nb [31] are used. The statistical un- 
certainty corresponds to the fit uncertainty on the area 
under the two Gaussian distributions, which is a fit pa- 
rameter. A systematic uncertainty of 3% covers the sys- 
tematic effects related to the luminosity and efficiencies. 
Contributions from background that peak at the J/tp 
mass are negligible. Using the electronic width Y J J^ = 
5.55 ± 0.14 keV, we determine the J/tp — > 7r + 7r~7r + 7r~ 
branching fraction to be: 



^J/lp^7T + 7r 7T + 7T 

= (3.67 ± 0.16 stat ± 0.08 svst ± 0.09 ext ) x 10~ 3 . (8) 



16 




1 2 3 4 0.5 1 1.5 2 2.5 0.5 1 1.5 2 2.5 




1VL (GeV/c 2 ) M 27I (GeV/c 2 ) M- (GeV/c 2 ) 



FIG. 12. Left: invariant tt + it~tt ± mass distributions (4 entries per event) for different regions in for the data without 
background subtraction (points) and signal simulation (histogram); Middle: tv + tv~ mass distributions (4 entries per event); 
Right: iy + tt~ mass distribution with another tv + tt~ mass in the p(770)° mass region. 



The external uncertainty (denoted ext) is dominated by 
the uncertainty of Y J Je [31]. This measurement agrees 
with the current PDG [31] value (3.55 ± 0.23) x 1CT 3 . 

A clear ip(2S) — > tt + tt~ /i + pr peak in the data of 
Fig. 3(a) is visible. Because the selection efficiency for 
TT + Tr~/i + fi~ is the same as for 7r + 7r~7r + 7r~, we can ex- 
tract the branching fraction ^(2S)^j/^.7r+7r- with J/ip — ¥ 
from a simple fit to the 7r + 7r~7r + 7r~ mass distribu- 
tion. Figure 13(b) shows the invariant mass distribution 
under the 7r + 7r~7r + 7r~ hypothesis in the ^){2S) mass re- 



gion. 

The measured width of the ip(2S) (« 20MeV) is, as 
in the case of the J/0 , dominated by the track momen- 
tum resolution, the intrinsic width of the ip{2S) being 
T^(2S) = 317 keV [31]. The effect of using the ir mass 
hypothesis in the [i track fit is negligible. The peak is 
described by the sum of two Gaussian distributions and 
the non-resonant 7r + 7r~7r + 7r~ contribution with a linear 
function. The area under the peak is used to determine 
the J/V' 7r+7r ~ branching fraction according to 
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FIG. 13. (a) Invariant 4ir mass distribution for the data 
around the J/ip peak with a fit, which consists of a sum of 
two Gaussian and a linear distribution, describing the peak 
and the non-resonant ir + ir~-K + ir~ contribution, respectively, 
(b) Invariant -k + tt~ -k + tt~ mass distribution, assuming all 
four particles are pions, for the data around the ip(2S) peak 
with a fit, which consists of a sum of two Gaussian and a 
linear distribution, describing the peak and the non-resonant 
it + -k~-k + -k~ contribution, respectively. 
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the following equation: 



M2S) 



and 



"i/>(2S)-*J/il>n+n- ' "jft>-*n r 

N(ip(2S) -> 7r+vr-/i+^-) 

dC/dE ■ due 
(84.7 ± 2.2 stai ± 1.8 syst ) MeV/c 2 nb 



(9) 



V»(2S) 



6tt 2 • dC/dE • e M c • G 
(49.9±1.3 stot ±1.0 syst )eV, 



where the latter result uses M, 



V(2S) 



(10) 
3686.09 ± 

0.04MeV/c 2 [31] and the conversion constant C — 
3.8938 x 10 n MeV 2 nb [31]. The statistical uncertainty 
corresponds to the fit uncertainty on the area of the two 
Gaussian distributions and the systematic uncertainty of 
3% covers the systematic uncertainties in the luminosity 
and efficiencies. Contributions from peaking background 
are negligible. The%^ +)1 - = 0.0593±0.0006 branch- 
ing fraction is known with very high precision. Using 
F V>(2S) 



M 4 „ (GeV/c") 



FIG. 14. Invariant n tt n n mass distribution for the data 
in the invariant mass range 3.2GeV/c 2 < Ma„ < 6.0 GeV/c 2 
(a) and 3.75 GeV/c 2 < M i7T < 4.25 GeV/c 2 (inset, b). 



C. Scan for additional resonances 

Fig. 14(a) displays the distribution for the data 
in the high invariant mass region. No clear signal can 
be identified. There is a hint of structure just above 
4GeV/c 2 . The inset, Fig. 14(b), shows this feature in 
more detail. 



VIII. SUMMARY 

In this paper, we present a measurement of the 
7r + 7r~7r + 7r~ cross section at effective center-of-mass ener- 
gies below 4.5 GeV, using e + e~ events with ISR collected 
in the vicinity of the T(AS) resonance. We achieve over- 
all uncertainties of 2.4% in the peak region defined by 



(2.38 ± 0.04) keV [31], the 4>(2S) -> J/^ir+n- 1.1 GeV/c 2 < < 2.2 GeV/c 2 , 10.7% below 1.1 GeV/c 



branching fraction is determined to be: 



B 



0.354±0.009 s 



t ±0.007 syst ±0.007. 



ext i 

where the external uncertainty is dominated by the un- 



certainty for r 



V>(2S) 



[31]. The measurement is slightly 



higher than the PDG [31] value B^ns^Mn+Tr- = 
0.336 ± 0.005, but agrees within the uncertainties. Our 
result is comparable in precision to the individual results 
used to determine the PDG average and agrees well with 
the most recent CLEO measurement B^s^M = 
0.3504 ± 0.0007 ± 0.0077 [45]. 



5.5% above 2.8 GeV/c 2 and 8.5% above 4.0 GeV/c 2 . These 
cross section results are much more precise than the cor- 
responding ones based on energy scans. 

The resulting contribution of the cr(e + e~ — > 
7T + 7r~7r + 7r~) cross section to the anomalous magnetic 
moment of the muon a l ^ ad (e + e~ — > ir + ir~ ir + tt~) is eval- 
uated according to the method described in Ref. [44] in 
the CM energy region between 0.6 and 1.8 GeV: 



3 /i 



'{tt+tt-tt+tt-) = (13.64 ±0.03 stat ±0.36 syst ) x lO" 10 . 

The cross section shows evidence of resonant sub- 
structure, with preferred quasi-two-body production of 
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Oi(1260)7r. There is an indication of a / 2 (1270)p(770) 
contribution to the final state. A detailed understanding 
of the four-pion final state requires additional informa- 
tion from states such as ir + TT~ir°ir° . 

The ISR events allow a study of J/ip and ip{2S) pro- 
duction. We measure the product of decay branching 
fractions and the e + e~ width of the J/i[> with the best 
accuracy to date, with the results: 



& J ftp— tTt + Tt 7T + 7T ' ^ ee 

= (20.4±0.9 stat ±0.4 S j, st )eV (13) 
= (49.9±1.3 atot ±1.0 syat )eV. (14) 
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